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How and when India collided with Asia is crucial for global climate and continental dynamics. We present 
new palaeomagnetic data showing that the Xigaze forearc basin of southern Tibet was located at 24.2 ± 5.9 0 N 
during 54-57 Ma, providing a direct constraint on the position of the southernmost margin of Asia at this 
crucial stage. Our study suggests 1) the age and locus of the initial India- Asia collision are at —50 Ma and 
~24°N, respectively; 2) Tibet resisted India's northward push during the first —16 Ma of initial impact 
from the collision and experienced little latitudinal displacement; and 3) Sometime a little after 34 Ma, 
Greater India was consumed and thicker Indian Craton subsequentiy made contact with Asia, resulting in 
~6° northward drift of Asia. Our model has implications for the process by which the high proto-Tibetan 
plateau formed and for the two slowdowns of India's convergence rate with Asia. 

The uplift histories of high plateaus are very important for understanding global climatic changes and 
continental dynamics. In particular, the Tibetan Plateau, which is well known as the 'roof of the world, 
has attracted much attention because of its high elevation, large size, and striking flatness 1 . Tibet mainly 
consists of three tectonic terranes: from north to south the Songpan-Ganzi, Qiangtang and Lhasa terranes, which 
accreted successively to Eurasia since the Early Mesozoic 2 . Later, during Cenozoic time 1 , the India- Asia collision 
occurred, contributing to the uplift of the Tibetan Plateau and also deforming large parts of Asia 2 4 . The uplift of 
the Tibetan Plateau may have triggered the trend of decreasing atmospheric C0 2 concentration and major change 
in global climate in Eocene 5 . 

The Yarlung Zangbo suture in south central Tibet marks the contact between the Lhasa terrane (Asia plate) to 
the north and the Tethyan Himalaya terrane (Indian plate) to the south (Fig. la). It is generally regarded as the 
remnant of the Neo-Tethys Ocean and consists of four major tectonic units (from north to south): the 
Cretaceous-Tertiary Gangdese magmatic arc, the Xigaze forearc basin, an ophiolite belt and the accretionary 
wedge containing Permian-Eocene rocks 6 . The Xigaze forearc basin formed at the southernmost continental 
(Andean) margin of Eurasia 7 and probably extended along the entire southern margin of the Lhasa terrane, 
representing the southernmost leading edge of Asia. Documented relics of the forearc basin at the present day are 
located between the subduction zone and the Gangdese belt and have a length of —660 km and width of —20 km 
(Fig. la). It was shortened by about 35-65% after the collision 8 , which is a minimum estimate because part of the 
forearc has been underthrust below the Lhasa terrane. The basin fill consists of Cretaceous-Palaeogene sedi- 
mentary deposits that were contributed from sources to the north in the Lhasa terrane 6 8 until the terminal forearc 
stage, which is represented by the Late Cretaceous-Early Palaeogene shallow marine sediments 6 . 

Despite more than three decades of active research, the lack of consensus regarding the timing of the col- 
lision 2,914 and the tectonic processes by which the Tibetan Plateau was raised 1 hinders understanding of the 
dynamic processes governing the kinematic history and structural development of Asia. Published palaeomag- 
netic studies have inferred inconsistent palaeolatitudes for the southern margin of Asia, which is the main cause 
for the wide range of collision ages (20-70 Ma) 9,12,13,15-21 , as well as post-collisional shortening estimates, which 
range from a couple hundred kilometers to —2000 km 919 . Van Hingsbergen et al. (2012) recently presented a 
compilation of palaeomagnetic data from Lhasa and Tethyan Himalaya terranes 13 showing a two-stage India- Asia 
collision at —52 and 25-20 Ma. Although this compilation is a timely service, the new tectonic scenario has been 
hotly debated 22,23 . These models use geophysical and geologic records between Lhasa and the Tethyan Himalaya 
terranes to constrain the India-Asia collision and lack a more direct measure of the collision, such as the coeval 
palaeolatitudes of Asia's southernmost and India's northernmost margins. 
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Figure 1 | Sketches of geology and sampling location for this study, (a), Geological map of the Western Indus-Yarlung Zangbo (IYZ) suture zone 
showing the position of Xigaze forearc basin and palaeomagnetic sampling area for this study (yellow star). Insert: present location of tectonic terranes in 
Tibet and surrounding regions; (b), Updated geological map of the study area by us, showing our palaeomagnetic sampling sites (A-B section) are located 
at the southernmost margin of the Asian plate. 



In an effort to provide constraints to resolve these questions, we 
carried out a geological and palaeomagnetic study of the Xigaze 
forearc basin in the Cuojiangding area of Zhongba County 
(29.9°N, 84.3°E) of southern Tibet (Fig. 1). The sedimentary 
sequences of the Xigaze forearc basin are characterized by the 
flysch-dominated Xigaze Group and shallow marine Cuojiangding 
Group (Fig. lb). The strike of well bedded Cuojiangding Group strata 
is approximate W-E. The Cuojiangding Group and the Xigaze 
Group strata were deformed by kilometer-scale synclinal folds, 
which are intersected by east-west trending Gangdese thrust struc- 
tures. The Gangdese thrust fault cut and placed the Xigaze and 
Cuojiangding Groups in the hanging wall and transported them 
southward over the Late Eocene-Miocene conglomerates 14,24,25 . The 
timing of slip on the Gangdese thrust has been inferred between 30 
and 23 Ma 25 . The Late Eocene-Miocene conglomerates were not 
folded. Thus, folding must have occurred during the Eocene, enab- 
ling a palaeomagnetic fold test on the age of magnetization. 

The Cuojiangding Group consists of the Padana, Qubeiya, Quxia, 
and Gyalaze Formations (Fig. 2). In the area we studied, the Quxia 
Formation is —80 m thick and is dominated by conglomerates with 
sandstone intercalations. The Gyalaze Formation is >120 m thick 
and composed of gray to dark gray foraminiferal sandy limestone 



and lithic sandstone. Well-lithified volcanic tuff layers occur near the 
middle part of the Gyalaze Formation, affording radiometric dates to 
further constrain the age of formation (54.8 ±0.7 Ma, Dr. Xiumian 
Hu, personal communication). Geochronologic investigations using 
U-Pb zircons 25 and fossils 6,26 also show that the age of the Gyalaze 
Formations is 50-60 Ma (Fig. 2). 

We collected 130 oriented palaeomagnetic samples from the 
Palaeocene Quxia and Eocene Gyalaze Formations of the upper 
Cuojiangding Group using a gasoline-powered drill. The main litho- 
logy of our samples is gray sandstone, sandy limestone and tuff. The 
samples were distributed along the northern and southern limbs of a 
syncline (Fig. lb). Where possible, samples were oriented using both 
magnetic and sun compasses. 

Results 

Rock magnetic results. A clear decrease in the heating curves during 
the temperature versus susceptibility experiments shows that all 
samples tested have well-defined Curie temperatures between 
560°C and 580°C (Supplementary Fig. SI). The susceptibility 
during cooling displays the same Curie temperature and is only 
slightly higher than that during the heating. These characteristics 
show that the main magnetic mineral is magnetite and does not 
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Xigaze forearc basin: Upper Cuojiangding Group 
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alter significantly during heating. Squareness (M rs /M s ) and 
coercivity ratios (H cr /H c ) obtained in the magnetic hysteresis 
measurements show that the mean magnetic grain size lies in the 
pseudo-single domain field of a Day plot 27,28 (Supplementary Fig. S2). 

Demagnetizations and palaeomagnetic directions. According to 
their demagnetization behavior, we assigned 4 grades to all the 
samples from the upper Cuojiangding Group (Supplementary 
Table SI, Fig. 3, and Supplementary Fig. S3). Out of 130 samples, 
31 samples were defined as grade A (= 'very good to excellent') 
(Fig. 3a-f). The natural remanent magnetization (NRM) of grade 
A samples is composed of two directional components. A high 
temperature component (HTC) is isolated that decays straight 
toward the origin, usually in the 450-590°C interval, indicating 



that the main magnetic carrier is magnetite and that this HTC is 
the characteristic remanent magnetization (ChRM). Several 
samples treated by alternating field (AF) demagnetization also 
display ChRM directions consistent with those from thermal 
demagnetization (Fig. 3b). Thirty eight samples were defined as 
grade B ( = 'OK and useful') (Fig. 3g-i). Their demagnetization 
behavior exhibited two or three directional components, but less 
clearly. The HTC component was observed at temperatures of 
300-580°C and decayed with a distinct trend toward the origin, 
but not as smoothly as for grade A samples. Only several HTCs 
required heating above 590°C to isolate the ChRM (Fig. 3g). Thirty 
three samples were defined as grade C ( = 'seems right but doubtful') 
(Supplementary Fig. S3a,b). Some of them exhibit a reversed 
intermediate unblocking temperature component (ITC) between 
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Figure 3 | Representative orthogonal projections of grades A and B palaeomagnetic demagnetization data from upper Cuojiangding Group plotted in 
in-situ coordinate. Solid (open) circles are the horizontal (vertical) plane projections. Scales are in mA/m. Straight lines represent the least-square-best 
fitting directions of the ChRM. See text for description of demagnetization behavior. 



300 and 500°C that is a little shallower but agrees approximately with 
the mean HTC direction of grades A and B (Supplementary Fig. S3a). 
Samples with directions that can only be constrained by great circles 
(Supplementary Fig. S3b) also were labeled as grade C. Twenty eight 
samples were defined as grade D ( = 'unacceptable'). These samples 
exhibited unstable directions that were unsuitable for applying 
principle component analysis to determine an HTC component 



(Supplementary Fig. S3c). After all low-stability components of 
NRM have been removed, the remaining HTC components are 
only accepted as the ChRM if they exhibit linear trajectories of 
demagnetization towards the origin and a reasonably consistent 
tilt-corrected direction from sample to sample that is significantly 
different from the present-day field in in-situ coordinates. We only 
used A and B sample directions for the final results in this study. 
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Table 1 | ChRM mean directions of upper Cuojiangding Group 
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Note: N, number of samples used to calculate mean direction; D g , l g (D s , l 5 ), directions in geographic (stratigraphic) coordinates; k g /k s , best estimate of Fisher's (1 953) precision parameter for the mean 
direction indicating the degree of clustering in geographic/ stratigraphic coordinates; a95 g /c(95 s , the radius of circle around calculated mean that has 95% chance of true mean lies with it in geographic/ 
stratigraphic coordinates. 



All low temperature components (LTCs) of the upper 
Cuojiangding Group are north and down and distributed around 
the present-day field direction (Dec = 0.6°, Inc = 46.4°). They yield 
a mean LTC which is obviously a modern-field viscous overprint (Dg 
= 0.2°, Ig = 44.6°, kg = 15.8, ot 95 = 4.0°, n = 85, in geographic 
coordinates; Ds = 181.4°, Is = 74.4°, ks = 3.9, a 95 = 9.0°, n = 85, 
after tectonic correction, see Supplementary Fig. S4). We calculated 
the mean ChRM directions by grade as shown in Table 1. The very 
few ChRMs with normal polarity in tilt-corrected coordinates were 
inverted to reversed. Grades A and B samples have nearly the same 
mean direction, whereas the mean inclination of grade C samples is 
slightly shallower. The grade A mean direction has the largest pre- 
cision parameter (k), and k decreases with the decline of reliability 
from grade A to C. To obtain a more accurate direction, all the C and 
D grade samples were excluded from further analysis. Finally, we 
remove the seven normal polarity directions from the formation 
mean because we suspect that they were remagnetized, albeit before 
tilting, as discussed below. 

Reliability of the ChRM directions and site mean direction. We 

sampled the upper Cuojiangding Group moving stratigraphically 
from bottom to top and each individual sample direction in a 
different sedimentary layer defines a palaeomagnetic site. In the 
end we obtained 69 site directions in the upper Cuojiangding 
Group. They are scattered in geographic coordinates, the 
inclination especially differing (up or down) in the two opposing 
limbs of the fold, and clustering increases after applying the 
structural correction. The ratio of the Fisher (1953) precision 
parameter, k s /kg = 2.0 is larger than the critical value of 1.25-1.39 
for 95% confidence 29 , indicating that the ChRM was acquired prior to 
folding. Furthermore, a Direction-Correction fold test of Enkin 
(2003) 30 for the site mean directions also gives optimal untilting at 
94.9±14.3% and a positive test. In addition, using the test of 
McFadden (1990) 31 , the fold test results are positive at more than 
99% confidence level (Supplementary Table S2): in in-situ 
coordinates, the value (22.043) of summation of cosine angles 
(SCOS), which was defined in Equation 14 of McFadden (1990) 31 , 
exceed the critical value (13.662) for the 99% confidence level. So 
there is a highly significant correlation between the ChRM directions 
and tectonics, and it is clear that the ChRM was acquired before 
tilting. Conversely, in the unfolded coordinates, the SCOS value 
(0.178) is far less than the critical value at 99% confidence level so 
the hypothesis that the ChRM was acquired before folding is tenable. 
Finally, according to geological evidence, the fold occurred during 
the Early Eocene, and thus we interpret the ChRM to be a primary 
remanence. 

The directions of 55 sites from the north limb of the upper 
Cuojiangding Group were used to calculate the virtual geomagnetic 
pole (VGP) latitudes as a function of stratigraphic height in Fig. 4. 
The magnetic polarity sequence is dominated by reversed polarity 
with only four normal samples scattered among the reversed sites. 
The tuff bed shown near the middle of the stratigraphic column is 
dated at 54.8 ±0.7 Ma, placing it squarely within C24r of the 



geomagnetic polarity timescale (GPTS) 32 , an entirely reversed polar- 
ity-chron spanning the interval 53.8-56.6 Ma. Although it is possible 
that the few normal samples might indicate three or four undocu- 
mented normal polarity subchrons, we think it more likely that they 
are isolated instances of pre-tilting remagnetization in a later normal 
field (e.g., Chron C24n). The minimum sedimentation rate for the 
upper Cuojiangding Group assuming that it lies entirely within C24r 
is ~7 cm/ka, which is reasonable for shallow sea facies. Therefore, all 
seven normal samples were removed from the population in further 
analysis, even though retaining them has no significant influence on 
the mean direction for the section (Table 1). As a result, we obtain the 
site mean direction for ~54-57Ma (red bold italic in Table 1 and 
Fig. 5): D s = 168.1°, I s = -42.0°, k s = 7.4, a 95 = 7.1°, n = 62 sites 
(after tectonic correction). 

The number of samples is less than ideal for the Elongation/ 
Inclination (E/I) method of Tauxe and Kent 33 , so we do not make 
any such correction. However, our palaeomagnetic mean inclination 
(42°±7.1°) is consistent with that (40.3°±5.6°) for the coeval vol- 
canic data from the upper Palaeocene to lower Eocene (56-47 Ma) 
Linzizong Group 13,17 . Agreement of our inferred palaeolatitude with 
Lippert et al. (2011)'s results from volcanic rocks 13-34 also supports 
this conclusion. For these reasons we conclude that our results from 
the upper Cuojiangding Group are not affected by inclination 
shallowing. 

The palaeomagnetic pole calculated from these tilt-corrected sedi- 
ments is located at 78°N, 329°E. The palaeolatitude of the sampling 
area (and thus the southernmost margin of this part of Asia) inferred 
from our palaeomagnetic data is 24.2°±5.9°N at —54-57 Ma. 
Compared with today's position of the sampling area, the southern- 
most margin of Asia has moved 5.7 ± 5.9° northward since that time. 

Discussion 

The motion of India relative to stable Eurasia can be reconstructed 
from the history of seafloor spreading in the Indian and Atlantic 
oceans and the rifting in east Africa. Various tectonic reconstructions 
have assumed that the dramatic change of India's northward con- 
vergence velocity between 50 and 35 Ma (dropping from —18 to 
—5 cm/yr) represents the onset of the India-Asia collision 9,10-35 37 . 
A more convincing, direct way to find out when India collided with 
Asia would be to determine when coeval palaeolatitudes of Asia's 
southern and India's northern margins first began to overlap. The 
overlap of the palaeolatitude curves derived from the Eurasian and 
Indian apparent polar wander paths (APWPs) 38 can't reflect the 
correct collision age because it occurs at less than 10 Ma (Fig. 6a), 
which is much too young. The likely reason is that shortening hap- 
pened within the Asia and India plates, and these palaeolatitudes are 
based on the present day shape of the continents. Thus to construct 
palaeolatitude curves for the Indian and Asian margins it is import- 
ant to use autochthonous palaeomagnetic data, that is data from as 
close to those margins as possible. 

Based on new compilations 12-13 , we review the palaeomagnetic data 
derived directly from the margin of Indian and Asian plates. Under 
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Figure 4 | Lithologic and magnetostratigraphic results from upper Cuojiangding Group in southern Tibet with declinations and VGP latitudes plotted 
as a function of stratigraphic level and showing the correlation with the GPTS 32 . 



the premise of meeting the foremost rigorous palaeomagnetic criteria 
(a proven primary signal and well-determined age constraint) 39 , we 
then add other palaeomagnetic reliability criteria (the correction of 
sedimentary inclination shallowing and averaging of palaeosecular 
variation (PSV) in volcanic-based data) to assess these individual 
poles. Together with our newest results, we extract robust palaeo- 
magnetic data to construct palaeolatitude curves of the Indian and 
Asian margins. 

On the Indian side, Late Cretaceous to Eocene palaeomagnetic data 
directly from the northern margin of India are scarce. Supplementary 
Table S3 lists available palaeomagnetic results from the Late 
Cretaceous Zhongshan and Palaeogene Zongpu Formations in the 
Dingri and Gamba areas of Tethyan Himalayas 40 43 . Cretaceous to 
Eocene marine strata in Dingri and Gamba areas of the Tethyan 
Himalaya represent the sedimentary development of the northern 



continental margin of the Indian plate 24 . Serious remagnetization in 
the Dingri area resulted in few sampling sites preserving primary 
remanence 12,17 ' 43 . Fortunately, recent palaeomagnetic results in the 
Gamba area of the Tethyan Himalaya (28.3°N, 88.5°E) yielded reli- 
able concordant palaeolatitudes for the Indian-northernmost mar- 
gin 42,43 . Patzelt et al. (1996) reported that India possessed equatorial 
palaeolatitudes of 4.7°S±4.4° (65-71 Ma), 7.5°N±4.8° (63-66 Ma), 
and 4.8°N±3.8° (55-63 Ma), respectively 42 , which were not affected 
by significant inclination shallowing 17 . Yi et al. (201 1) obtained higher 
resolution results from the same section with palaeolatitudes ranging 
from 6.6°N±3.5° (59-62 Ma) to 11.1 0 N±2.5° (56-59 Ma) in the 
lower and upper Zongpu formations 43 . Positive fold and reversal tests 
indicated that the ChRM directions are primary remanence. 

On the Asian side, Late Cretaceous data stem from studies of the 
Takena and Shexing Formations in southern Tibet (Supplementary 
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Figure 5 | Equal area stereonet projections of palaeomagnetic site directions from the upper Cuojiangding Group with a site mean direction (red black 
star, 95% confidence interval indicated in blue shadow): (a), in-situ site directions; (b), tilt-corrected site directions. Solid black circles indicate lower 
hemisphere projections and open black circles are directions in the upper hemisphere. 



Table S4) 15,19,21,44 - 46 , with palaeolatitudes ranging from ~12°N to 
~25°N 12 . Basically, all these studies more or less agree with Tan 
et at (20 10) 19 if large inclination shallowing is assumed. Tan et at 
(2010) studied the Late Cretaceous redbeds of the Shexing formation 
(65-110 Ma) and volcanics intercalated within them (85-110 Ma)in 
the Linzhou area of southern Tibet (29.9°N, 91.2°E). The ChRM 
directions of both the redbeds and volcanics pass the fold test at 
99% confidence and can be interpreted as primary remanence. 
Importantly, this study is the only one with enough intercalated lava 
flow samples to do a reasonable job of averaging PSV. The calculated 
pole dispersion Svgp = 14.9° is consistent with known dispersion 
models for PSV: S = 14.5° in Johnson et at (2008) 47 and S = 12.5° in 
Biggin et at (2008) 48 for — 21°N of palaeolatitude, which is in agree- 
ment with that for the coeval redbeds corrected by the E/I method 33 . 
This proves that E/I analysis correction is needed and works well for 
the Late Cretaceous redbeds of the Shexing formation. Other studies 
agree with Tan et at (2010) when the inclination is corrected by E/I 
analysis 13 . To sum up, a palaeolatitude for the southern Asian margin 
of — 25°N is remarkably consistent for the reference point (29.9°N, 
84.3°E) during Late Cretaceous time. 

Most of the Tertiary results are from the Palaeogene Linzizong 
volcanic rocks in the same area on the Lhasa terrane (Supplementary 
Table S4). The Linzizong Group was classified into three formations 
(Dianzhong, Nianbo and Pana Formations) with accurate geo- 
chronological constraint ranging in age from 38 to 65 Ma 49 . 
Palaeolatitudes derived from Linzizong Group are highly inconsist- 
ent and controversial 12,13,15-20,45 . Nevertheless, Dupont-Nivet et at 
(20 10) 17 and van Hinsbergen etat (20 12) 13 compiled and synthesized 
palaeomagnetic results from the Linzizong volcanics, making the 
case that the results probably averaged out the PSV because of the 
large numbers of volcanic sites, and obtaining an average palaeola- 
titude quite close to that from our sediments. 

Palaeomagnetic data from numerous Late Eocene lava flows 
(—34 Ma) of the Qiangtang terrane 34 also appear to have averaged 
out PSV reasonably well, the calculated pole dispersion agreeing with 
that expected for palaeosecular variation 47 . Their palaeomagnetic 
pole gives a palaeolatitude of 28.7±3.7°N for a reference point at 
— 33°N, ~88°E. When transferred to the southern margin of Asia 
(29.9°N, 84.3°E), these Eocene volcanic results suggest a palaeolati- 
tude of 25.6±3.8°N for 32-38 Ma assuming that the amount of N-S 
crustal shortening or rotation within the Lhasa terrane during 
Neogene time can be ignored 13 . Compared with our results, the 
palaeolatitude difference of Tibet between —50 Ma and —34 Ma is 
1.4±5.9°, which indicates negligible to at most modest latitudinal 



displacement of southernmost of Asia considering the error. Thus, 
this synthesis of recent palaeomagnetic results including our own 
data described above, derived from different parts of Tibet and 
depicted by the red dashed line in Fig. 6a, shows that the palaeola- 
titude of the southern margin of Asia was remarkably consistent 
around 24-25°N, moving almost imperceptibly northward from 
Late Cretaceous to Late Eocene time. 

Across the suture zone, palaeolatitudes for the northernmost 
Indian margin increase progressively with decreasing age from 
4.3°S±4.4° at -65-71 Ma, to 7.2°N±3.5° at -59-62 Ma, to 
11.9°N±2.5° at -56-59 Ma for the reference point (29.9°N, 
84.3°E) 12,13,17,42,43 . This northward drift of the Indian margin, indi- 
cated by the green dashed line in Fig. 6a, is not exactly parallel to the 
blue solid line derived from the Indian APWP (Fig. 6a). If significant, 
the divergence in the lines would imply that extension was occurring 
in continental India during the Cenozoic until the time of collision. 
Van Hinsbergen et at (2012) champion the role of very large exten- 
sion commencing in the Early Cretaceous 13 . We do not choose to 
speculate beyond the Cenozoic time frame of our study, but simply 
note that the discrepancy in the two palaeolatitude lines indicates 
that the margin of Greater India at the time of collision lay 
— 1500 km or more to the north of the Indian Craton than it does 
today. 

In Fig. 6a extrapolation of the palaeolatitude drift green-dashed 
line shows that the Indian margin at the 56-59 Ma time mark still 
needed to move —13° farther north before reaching the Asian mar- 
gin (indicated by the red dashed line). The intersection of the two 
palaeomagnetically determined drift lines gives the time of initial 
collision as well as its palaeolatitude. Using this configuration for 
pre-collisional Asia and India and the uncertainties shown in 
Fig. 6a, the best estimate for the age of initial collision is 49.7 Ma 
with upper and lower 95% confidence bounds of 44.3 and 55 Ma. 
The Gyalaze Formation (Fig. 2) marks the cessation of the marine 
strata along the suture zone and has been dated at —50 Ma 6 , suggest- 
ing that the ocean between India and Tibet was closed and contin- 
ental collision had started. Therefore, the —50 Ma onset collision age 
is constrained by both palaeomagnetic and sedimentary studies from 
the same strata. 

Based on this kinematic and other supporting evidence, we pro- 
pose a tectonic model for the India- Asia collision (Fig. 7). The north- 
ern edge of Greater India was around latitude 7°N at —60-62 Ma 
(Fig. 7a) 43 . The widespread Cretaceous limestone and sandstone 
covering most of southern Tibet 1 indicate that it was a stable platform 
when these sediments were deposited. At around 57.5 Ma, Greater 
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India was drifting northward at a velocity of —18-19.5 cm/yr 10 and 
the palaeolatitude difference between the Indian and Asian margins 
was —13° (Fig. 7b). At —50 Ma, the palaeolatitude difference 
reached zero at — 24 °N, marking the initial continent-continent col- 
lision between India and Asia (Fig. 7c). With estimated subduction 
rates of —7-8 cm/yr 35,37 , the thinner non-cratonic part of Greater 
Indian margin (with a width —1500 km) could have been almost 
entirely consumed by —34 Ma. Thickened and uplifted Tibet 1 
resisted thinned Greater India's northward subduction 50 during the 
first —16 Ma of initial impact from the collision and experienced 
little latitudinal displacement, as reflected by the palaeolatitude data 



(Fig. 7d). Probably by the end of the Eocene (—34 Ma) or a little later, 
the thicker Indian Craton 50 and Asia made contact, increasing com- 
pressional stress on the Asian plate 50 and resulting in —6° shortening 
or northward drifting of Asia to the present position (Fig. 7e). This 
model better explains the onset of Indochina extrusion by —20- 
35 Ma 51 and may also have led to the two reductions in convergence 
rate of India with Asia at 45-40 Ma and at 20-10 Ma deduced by 
Molnar and Stock (2009) 37 . 

Our model agrees in age and palaeolatitude of collision with van 
Hinsbergen et al.'s model for Tibetan-Himalayan microcontinent- 
Asia collision 13 , but differs considerably in the pre-collisional 
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Figure 7 | Model showing the kinematic and dynamic evolution of the India- Asia collision. See text for details. 



configuration of Greater India and post-collisional kinematic recon- 
structions. We describe these differences below. 

Because no suture zone has been identified to suggest Tethyan 
Himalaya was ever a terrane separate from Indian continent 52 , recon- 
structions 42,43 based on the same Palaeocene results from the Tethyan 
Himalaya used in our model have proposed a large northern extent of 
continental Greater India (—1500 km) since Late Cretaceous. In 
addition, new palaeomagnetic results from the Tethyan Himalaya 
also suggest >1300 km northern dimension of Greater India since 
Triassic 52 . Thus, in our model we start with a large continental 
Greater India for the pre-collision configuration of India. In contrast, 
by analyzing the Early Ordovician to Palaeocene palaeomagnetic 
poles for India, van Hinsbergen et al. (2012) 13 propose that 
—2500 km of N-S Cretaceous extension occurred within Greater 
India while the Tibetan Himalayan microcontinent rifted away from 
India in the Cretaceous by opening an oceanic Greater India Basin. 
Furthermore, in our model the palaeolatitude difference of Tibet 
between —50 Ma and —34 Ma is 1.4±5.9°. This probable near- 
standstill in palaeolatitude suggests that the continental convergence 
was mainly accommodated by the subduction and consumption of 
the Indian continental crust during this interval. Van Hinsbergen 
et al, on the other hand, propose in their model that subduction of 
their extended oceanic Greater Indian Basin continued significantly 
later, until -23 Ma. 

In our kinematic reconstructions, the convergence history can be 
matched by numerical models 50 showing that two stage subduction 
of the denser thinned Greater India and thicker cratonic India and 
are in agreement with geological and geophysical observations. 



India's northward convergence as indicated by analysis of sea-floor 
magnetic anomalies slowed markedly around Anomaly 21 -22 9,10,35,37 
(48-50 Ma, Fig. 6b-A). The Eocene age of the last sedimentary mar- 
ine facies dated as foraminiferal zones P8-P12 across the Himalayan 
belt 11 (Fig. 6b-B) correlates well with the initial stage of collision. 
Evidence for subduction-related magmatism 1 (Fig. 6b-C) and the 
first arrival of sedimentary material from Asia onto the Indian con- 
tinental margin 53 (Fig. 6b-D) also corresponds with the timing of 
initial India- Asia collision at —50 Ma. In addition, ultrahigh-pres- 
sure rocks formation and exhumation currently exposed in the 
Himalayas dated between 57-50 Ma 54 (Fig. 6b-E), suggesting the 
onset of Indian continental subduction below Asia during this time 
interval. 

The standstill feature of the Tibetan Plateau in our model during 
the first — 16 Ma of the collision is also consistent with some geologic 
and geophysical evidence in the region. During —50-34 Ma, the 
existence of thickened high proto-Tibetan plateau 1 , the significant 
eastward shift (—40°) of the Euler pole for India, and the rearrange- 
ment of the plate boundary in the Indian Ocean 35 all favor the hypo- 
thesis that a large part of Greater India thrust beneath Tibet and 
thickened the crust by mechanical coupling 55 . 

In conclusion, the kinematic model, based on our palaeomagnetic 
results establishing the palaeolatitude of the southern margin of Asia 
and on other palaeomagnetic and geological evidence, shows that 
Greater India was converging with the southern margin of Eurasia 
in Tibet at a palaeolatitude of — 24°N and began to collide at around 
50 Ma. For about the next 16 Ma the southern margin of Asia 
remained almost fixed while thinner Greater India was subducting 
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and causing shortening in the Himalaya, early uplift of the Tibetan 
Plateau, and reorganization of plate boundaries in the Indian Ocean. 
Sometime a little after 34 Ma, Greater India was consumed and 
thicker Indian Craton subsequently made contact with Asia, result- 
ing in —6° shortening or drifting northward of Asia. The model also 
provides a link between the timing of uplift of the Tibetan Plateau 
and global palaeoclimatic changes, namely for the Cenozoic global 
cooling that started at —50 Ma 5 and persisted during the Eocene- 
Oligocene transition (—34 Ma), as mountain building process and 
physical weathering of high plateau consume large quantity of C0 2 5 . 

Methods 

All demagnetization treatments and remanence measurements were performed in a 
shielded room with residual fields less than 300 nT at South China Sea Institute of 
Oceanology. Most of the samples were treated with stepwise thermal demagnetiza- 
tion. Some samples were also progressively demagnetized using alternating fields. The 
remanence after each step was measured on a 2G cryogenic magnetometer. Rock 
magnetic experiments were performed at the University of California, Santa Cruz. 
Magnetic hysteresis curves were measured on a Princeton Measurements MicroMag 
2900 alternating gradient. The temperature versus susceptibility curves were mea- 
sured using a Kappabridge KLY-2 susceptibility bridge with a CS-2 furnace. 

Palaeomagnetic directions were identified using principal component analysis to 
obtain a best-fitting demagnetization trajectory 56 . A few samples were analyzed with 
great circle methods 57 . The formation means of best-fit directions were calculated 
using Fisher statistics 58 , while those that included some great circles were calculated 
following McFadden & McElhinny (1988) 59 . Confidence limits were calculated 
according to Coe et al (1985) 60 . 
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